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Abstract 

By using slave particle (slave boson and slave fermion) technique on the Bose-Hubbard model, 
we study the finite temperature properties of ultracold Bose gases in optical lattices. The phase 
diagrams at finite temperature are depicted by including different types of slave particles and the 
effect of the finite types of slave particles is estimated. The superfluid density is evaluated using 
the Landau second order phase transition theory. The atom density, excitation spectrum and 
dispersion curve are also computed at various temperatures, and how the Mott-insulator evolves 
as the temperature increases is demonstrated. For most quantities to be calculated, we find that 
there are no qualitatively differences in using the slave boson or the slave fermion approaches. 
However, when studying the stability of the mean field state, we find that in contrast to the slave 
fermion approach, the slave boson mean field state is not stable. Although the slave boson mean 
field theory gives a qualitatively correct phase boundary, it corresponds to a local maximum of 
Landau free energy and can not describe the second order phase transition because the coefficient 
<i4 of the fourth order term is always negative in the free energy expansion. 

PACS numbers: 03.75.Lm,67.40.-w,39.25.+k 
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I. INTRODUCTION 



Ultracold atoms in optical lattices with highly tunable parameters have provided an 
unique opportunity to simulate strongly correlated phenomena in condensed matter physics 
DDD . The dynamics of such Bose atoms can be described by a Bose-Hubbard model, and 
it was predicted that there would be a quantum phase transition from the superffuid to the 
Mott-insulator phase induced by varying the depth of the optical potential ^J. Recently, 
this phase transition has been perfectly realized by Greiner et al. by means of loading 
87 Rb atoms into a three-dimensional optical lattices. From then on the Bose-Hubbard model 
to the cold atomic gas in optical lattices has received extensive studies. 

The Bose-Hubbard model, with on-site interaction and tunnelling between nearest neigh- 
bor sites, was US ed by FishOT e, * fl to investigate the boaoas ia periodic and/e, raadoaa 
external potentials. There exist two types of phases, the superffuid and Mott-insulating 
phases, in this model at an integer lattice filling fraction and zero temperature. In the su- 
perffuid phase, the tunnelling term dominates and all atoms occupy the identical extended 
Bloch state. The ground state of this system can be well described by a macroscopic wave 
function with long-range phase coherence. In the Mott-insulating phase, the interaction 
dominates and the ground state of the system instead consists of the localized atomic wave 
function without phase coherence ^|. This Mott-insulating phase can be characterized by 
an integer filling factor, the existence of a gap for particle-hole excitation, and zero com- 
pressibility. By varying the parameters such as the density and the external potential, the 
system would undergo a quantum phase transition and evolve from the superffuid phase 
to Mott-insulating phase. In the past, various theoretical approaches have been used to 



coupling expansion 



,0], Gutzwiller projection ansatz jj, 
Carlo simulations |!3l. Il4j. and other mean- field approximations [5 
there are less studies focusing on the nonzero temperature properties 
and this will be the main topic of this paper. 
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Slave boson technique was developed to deal with the interacting fermion systems for 
the convenient use of mean-field approximations j^]. It has been extended to study the 
hard-core bosons on a lattice in a functional integral representation by Ziegler j^j , in which 
there are only two states per site: empty and singly occupations. A similar formulation for 
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interacting boson systems has been obtained by Fresard 24]. The analytical expression for 
the Mott-insulating lobes was obtained and the density- density correlation function was cal- 
culated. Recently, this approach has been adopted by Dicherscheid et al. jl?! to investigate 
the nonzero temperature behaviors of the ultracold atoms in optical lattices. These authors 
showed that the calculated phase diagram at zero temperature agrees well with the results 
by other mean field methods |(| and they further depicted the phase diagrams at finite 
temperatures. The slave boson and slave fermion approaches are equivalent in principle 
without any approximation. However, the difference comes into view if some mean field ap- 
proximations are made in practical calculations. Bearing this in mind, one of us with Chui 



21( performed a comparative study between the slave fermion and slave boson approaches 



under the same mean-field approximation. It was found that both approaches give the same 
qualitative phase diagram, but the quantitative behaviors by slave fermion approach are 
more accurate. This may be related to the Fermi statistics which automatically excludes 
two same types of slave fermions from occupying the same site when the constraint, one site 
can be occupied by only one slave particle, is relaxed in the mean field approximation. 

In this paper, we will extend the slave particle technique to investigate the various finite 
temperature properties of the Bose-Hubbard model. We will study the finite temperature 
phase diagram, atom density, superfluid density, and excitation spectrum for various tem- 
peratures and system parameters. The comparison between slave fermion and slave boson 
approaches will be performed. We show that although the mean field theory of slave boson 
gives a qualitatively reasonable phase boundary, the negative coefficient 04 of fourth order 
term in the Landau free energy expansion means that the zero of the derivative of the free 
energy is a local maximum and thus the slave boson mean field state is not stable, i.e., the 
slave boson mean field theory is unable to describe the second order phase transition of the 
normal liquid/superfluid. This shortcoming of the slave boson approach does not exist in 
the slave fermion approach. In this sense, using the slave fermion approach to study the 
strongly correlated Bose atoms may be preponderant. 

This paper was organized as follows. In Sec. |H] we review the basic formalism of the 
slave particle technique to the Bose-Hubbard model. In Sec. IIHI we focus on the critical 
temperature of superfluid-normal phase transition, and depict the phase diagram. The atom 
density and the compressibility are calculated in Sec. IIVI In Sec. El the superfluid density 
is determined and the stability of the mean field states of the slave particles are discussed. 
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Sec. |VT] is devoted to calculate the excitation spectrum. Conclusions are given in Sec. I VI II 
In Appendix A, the details of the calculation of the coefficient 0,4 is provided. 

II. SLAVE PARTICLE FORMALISM OF BOSE-HUBBARD MODEL 
A. Slave particle technique for a bosonic system 

For a bosonic system, the creation operator a\ and annihilation operator on site i can 
be defined as follows in the occupation-number representation: 

a\\a)i = Va + l\a + 1);, 

ai\a)i = y/a\a - (1) 

They obey the basic boson commutation relation: [ai,aj] = 5{j. The state \a)i is an eigen- 
state of the particle number operator Ni = a\ai, which counts the number of bosons on site 
i, with the eigenvalue a. On a single lattice site, the occupation number can be any non- 
negative integer. Thus the boson creation and annihilation operators can be decomposed 
into 

a\ = \l) u (0\ + V2\2) u (l\ + .-. = J2 Va + l\a + l)«(a|, 

a=0 

a< = |0)«(1| + V2|l)«(2| +-. = X; V^+l|a>«<Q! + 1|- (2) 

The justification of this decomposition is that it does satisfy the original commutation 
relation [aj,aj] = 8ij. Next, we identify every occupation state on a site as a type of slave 
particle, that is to say, mapping and i(a\ to the creation operator a ai and annihilation 
operator a a ^ of the slave particle. In a slave fermion approach these operators are forced 
to satisfy the fermion anticommutation relation {a a ,i, a ^j} = dapfiij, and in a slave boson 
approach they are forced to obey boson commutation relation [a a ,i,a^j] = &ap&ij- Then a\ 
and di can be rewritten as 

at = Va + lal +l i a a>h 

a=0 (3) 
a=0 
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When inserting the above equations in [cij,at], one finds that, only when the constraint 

^2 n i=^2 a U a a ,i = h ( 4 ) 

is satisfied, could the original boson commutation relation [a^aj] = 5^ be reproduced by 
either slave fermion or slave boson approaches. This implies that the slave particle trans- 
formation Q along with the constraint (jlj) and {a a> i,a^j} = 5 a pbij (slave fermion) or 
[ a a,i, a ^j] = 8 a p8ij (slave boson) can describe the whole physics of the original boson sys- 
tem. 



B. Bose-Hubbard model and functional integral representation 

In the second quantization, the translationally invariant many-body Hamiltonian of cold 



3 , 



Bose gases confined by an external optical lattice potential V ext {r) is given by 

H = f dr¥(r) f-|^V 2 + V ext (r) - fi) *(r) 

+ - / drdr'^ t (r)^ t (r')l / (r-r , )^(r')^(r), 

where ^(r) (^'(r)) is the boson field operator that annihilates (creates) a particle at the 
position r, V(r— r') is the two-body interatomic potential, and /z is the chemical potential. In 
the case of a dilute cold atom gas, we can approximate V(r — r') with an effective interaction 
g5(r — r'), where g = 47ra s h 2 /m with a s the s-wave scattering length and m the mass of 
the atoms. When expanding the field operators in the Wannier basis and keeping only the 
lowest vibrational states, namely ^(r) = ^2 i a i w(r — rj), eq. (J3J) can be rewritten as the 
Bose-Hubbard Hamiltonian Q]: 

H = -t^alaj- ft^2n i + -^2ni{n i -i), (6) 

<ij> i i 

in which rii = a\ai is the particle number operator. The symbol (ij) denotes the sum 
over all nearest neighbor sites, t and U are the hopping amplitude and on-site interaction, 
respectively, 

t = J drw*(r - n) ("|^V 2 + V ext (T)^j w(r - r,), 
U = g [ dr\w(r)\ 4 . (7) 



For small occupation per site, one can use the single particle Wannier function to calculate 
these parameters. However, in the case of the multi-occupation, these parameters have 
to be calculated by considering the interaction broadening of the Wannier function 27]. 
Substituting eq. (j2J) into eq. (JHJ), the Bose-Hubbard Hamiltonian can be replaced by 



H = ~ t + + l4 +lii a Q)i aJ ^ap+ij 

<ij> a, (3 

u 



(8) 



Following the steps in Refs. |28|, |29|], we write the partition function as an imaginary time 
coherent state path integral: 

Z = Tre~P H = I Da n Da n DXe- s[aa ' aa ' x] 



f 13 { U 
S[a a ,a a ,\]= / dr\ a a ,i[d T — afi + — a(a — 1) — i\i]a c 

+i y \ - 1 Va + lv 7 /^ + la<a+l,i a a,iQ>l3,ja'f3+l,j (, 

i <ij> a,f3 



(9) 



where a a j and Cl a i £1X6 introduced as ordinary complex numbers in the slave boson approach, 



and as Grassmann variables in the slave fermion approach satisfying t 
i.e., {a a ,a a } = {a a ,a a } = {a QJ a a }, J da a l = and J da a a a = 1 
integrals over them are 



re Grassmann algebra, 
The gaussian 



/ Da a Da a exp < — a a A a pap > = (detA)^ 1 = e 

J I a/3 J 



±Tr [In A] 



(10) 



where ± correspond to integral over Grassmann variables or ordinary complex numbers. 
The Lagrange multiplier field A,(r) comes from the constraint (@J, namely, FJ i $(Yla n i ~~ !)• 
The unit has been set to h = kg = 1 in all formulas. In order to decouple the hopping term, 
a Hubbard- Stratonovich transformation is performed by adding a complete square term to 
the action which contributes to the partition function a constant, 



exp 



<ij> a 

~ Y V« + la aJ a a+1 j) 



(11) 
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Then the partition function is replaced by 



D$*D<&Da n Da n D\e- s eff^^ x \ 



S e ff[$, a a , A] 



*■ EE' 



d T — a\i + — a(a — 1) — iX L 



+i x i + * J2 (®i $ i ~ ®i V« + la a>j a a+ljj 

i <ij> a 

a ) 



;i2) 



The Hubbard-Stratonovich field <3?j introduced here can be identified as the order parameter 
of superfluid for (<3>j) = (J2 a V a + lda,i a a+i,i)) = We then perform a Fourier transform 
on the field Ac 

i 

3 i(k-ri— cj„t) 



-4; 



^E-4 



(13) 



k,n 



where L is the total number of sites in the optical lattice, u n is the Matsubara frequency, 
which equals (2n + l)ir//3 and 2mr/f5 for fermion and boson fields, respectively. The effective 
action S e ff now reads 



k,n o 



/ 

-iuj n — an + — a[a — 1J 



^a,kn 



k,n 



k,q, n,n' a 

— ^2 ^2 FTa { ^ a ^ l^a+l,(k+k')(n+n') aQ . k «) ^k',n' 



(14) 



k,k',n,n' a 



+ (y/a + la a ,kna a +l,(k+k')(n+n')) ^k',n' 

where et = 2£ Y^ =1 cos(fcja) with d and a being the dimension and spacing constant of the 
lattice. So far, we have obtained an effective action in a functional integral representation 
by the slave particle approach. It is a reexpression of the original Bose-Hubbard model since 
all transformations we made are rigorous. This effective action is the starting point of our 
calculations. 



C. Perturbation theory 



The system can not be exactly solved with the constraint (@J), which means there is 



exactly one type of slave particle per site (22, 



301 ] . We then relax it to one slave particle per 



site on average over the whole lattice. To realize this, one can replace all A ki „ in the action 
(fHj) with a constant An n - This approximation is widely used in dealing with the interacting 
fermion systems 0, [icj • It implies that the multi-occupation of the slave particles on one 
site is allowed. However, the behaviors of the slave boson and slave fermion are slightly 
different in the relaxing process. The multi-occupation of the same type of slave boson is 
allowed, while the multi-occupation of the same type of slave fermion is forbidden by the 



2l\ . This will lead to some quantitative differences in the results by these 



Pauli principle 

two approaches. 

According to the Landau phase transition theory, the order parameter $ near the critical 
point is small and then the perturbation theory can be used. In the following, we will try to 
integrate the slave particle field out of the action (|14|) . and perform perturbation calculation 

After relaxing the constraint (@J), the action (fTlj) can 



toward $ along the way in Ref. 
be divided into two parts: 



S e ff[$,a a ,\] = S + Sj 
S = iL(3\ + y^y^ a a ^ n [-iu n + c(a)] a aM + ^ e k |$ kin | 2 = Sq P + e k |$ k)ri 

k,n a k,ro k,n 

Si = — Pf-p. { + ^a+l,(k+k')(n+n') a a,kn J $ k ', n ' 



k,k, n,n' a 

+ (Va + Ia ajkn a a+li ( k+k ')( n+n /)^ 3> k ' n /j , 
in which A = An,o/\/I/3 is a constant and 

c(a) = —iX — apt + a (a — 1)77- (15) 

The partition function Zq of non-interacting slave particles comes from the contribution of 
the zeroth-order term Sq P and is given by 

Z = e~ pn ° = [ Da a Da a e- s ° P (16) 



where Qq is the zeroth-order thermodynamic potential and has the form: 

- Q = ± - ^ln(l ± e -*W) f (17) 

a 

where the +(— ) sign corresponds to the slave fermion (slave boson), respectively. Then we 
can define the average of an operator A with respect to Sq P as 

(A) = Y Q J D aaDa a A{a a ,a a )e- s ° P . (18) 



For small $, e efi can be expanded in terms of Si [17, 29], i.e., 



(19) 



After integrating out the slave particle field a a and a a , we arrive at a new effective action 

S E,eff, 

1 



l-(5/)o + 7TW)( 



(20) 



It is straightforward to calculate (Sj)o and (S])o using the Wick's theorem and gaussian 
integral formula (eq. (fTUJl ). 



(S/) = 0, 



(21) 



k,k', n,n' 



X a r 



a+l,(k+k')(n+n') tt a+l,(k+k')(n+n') 



>« 



£ 



k,k', n,n' 



f|^kV 

L/3 



-za;„ + c(a) 



x- 



(22) 



-iuj n+n , + c(a + 1) ' 

where =F correspond to the slave fermion and slave boson. After performing the sums over 
Matsubara frequency u n and k, we reduce {S]) to 

<^>0=2E4l$k,n| 2 $> + l)- 
k.n a 

where n a is the occupation number and equal to 

1 

n c 



n a - n a+1 



-iuj n — /i + all' 



(23) 



(24) 



exp{(3[-i\ -an + a(a - l)U/2]} ± 1 ' 

in which the + and — sign correspond to slave fermion and slave boson respectively. Because 
(Sj) is small, we have 



l-(S/>o + ^/>o 



(25) 



and can rewrite eq. (|2U|) as 



e SE,eff _ e -/3^0-Ek,„^k|-fk,n| 2 + ^S?>0_ 



(26) 
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Then the effective action of the order parameter field SE,eff[<&*, ma Y be obtained 

S E)tiff [&,$\ = PQ -J^K^iKi^Kn, (27) 

k,n 

where the Green's function G(\c,iu n ) is defined by 

-G-\k,iu; n ) = e k - e* V(a + 1) """"^ .. . (28) 
which will be the starting point of our analysis. 

D. Mean-field approximation 

In terms of eq. (|27j). the partition function can be calculated by gaussian integral 

Z =e -/3n = y D$D$*e~ SjJ ' e « 

D$D$*e^° 0+Ek >" *£, n G- 1 (k,t«„)* k , B (29) 

=e -/3n -E k ,„ln[-/3G- 1 (k,io; n )] ) 

from which the thermodynamic potential Q can be extracted: 

Q = + \ 2 ^h^G-^k, iw n )]. (30) 

^ k,n 

We then perform a saddle point approximation to the constraint field A, which means that 
we only choose the A minimizing the thermodynamic potential: dQ/dX = 0. In addition, 
the particle number conservation condition requires —dQ/d^i = N. When inserting eq. (|30|) 
into these two conditions, we have 

L^_J\-n ) - - 2^ G(k, iu n ) — = 0, 

(31) 



L^an +-^G{k,iu n ) — = N. 



u k,n 

The mean-field approximation means the last terms of the above equations may be neglected. 
That is, all the fluctuations coming from the Green's function would not be considered. Then 
the following two equations can be derived, 



5>° = 1, (32) 



a=0 
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^cm Q = — = n, (33) 

where the n = N/L is the average particle density. One can see that eq. f!32[) . which implies 
there is one slave fermion per site on average, is exactly the relaxed constrain (J3J. 

III. THE SUPERFLUID-NORMAL PHASE TRANSITION 

As we know, the quantum phase transition from superfTuid to Mott-insulator only occurs 
at T = 0. At finite temperature, a superfluid-normal phase transition will be induced 



17. 18. 



21 



231 ] . The Landau theory shows that, near the critical point of superfluid-normal 
phase transition, the order parameter of superfluid $ is small and the Landau free energy 
can be expanded in terms of it. The critical point can be determined by the coefficient of 
the second order term |$| 2 , which is Cr _1 (0, 0) = in our case According to eq. 



28]) and noting that e = zt, we have 



n a+l 



o + l- -=- = 1, (34) 

where p, = fi/zt and U = U j zt are dimensionless chemical potential and on-site repulsion 
strength with z being the number of nearest neighbors. The dimensionless critical temper- 
ature T c = T c /zt can be obtained by solving eqs. (J2H), (JH2J, (]3*3^1. and (jHjj) together. Note 
that our results have a clear mean field nature and the dimensionality of the lattice appears 
only as a numerical factor. Because there are infinite types of slave particles, a cut-off toward 
a should be made when solving these equations. It has been shown by one of us and Chui 
JQ that, for a large U a small cut-off (e.g. q>m = 3) is a good approximation, but in the 
region of small U the contributions from large a slave particles cannot be neglected and a 
larger is required to obtain a quantitatively reliable result. In this work, we will work at 
a relative large olm ( say, 9) to obtain a more accurate phase diagram of superfluid-normal 
transition and estimate the effect of the finite types of slave particles. 

When considering the commensurate state and restricting the average density to n — 
N/L = 1, we obtain the superfluid-normal phase diagrams with various cut-off «m in both 
slave fermion and slave boson approaches. It can be seen from Fig. Q that, in both ap- 
proaches, the curves of various cut-off are different from each other in the small U region. 
As U increases, these differences become small and disappear gradually when approaching 

11 




FIG. 1: The superfluid-normal phase diagrams with various cut-off ot\j in the slave fermion ((a)) 
and the slave boson ((b)) approaches. The average density is set to be n = 1. T c = T c /zt and 
U = U/zt are the dimensionless critical temperature and on-site interaction, respectively. 



the critical point of superfluid-Mott insulator transition. This result gives a further sup- 
port to the statement of Ref. 2l| that more types of slave particles should be taken into 
account in the small U region. Moreover, we find that the T c — U curves with olm > 7, 
which are not shown in Fig. ^ almost coincide with the curve with au = 6. This implies 
that even in the small U region the contributions from the slave particles with a > 6 can be 
reasonably neglected. To make the effect of finite cut-off au clearer, we plot the position 
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FIG. 2: The position of local maximum U max (panel (a)) and the critical temperature T c $ at U = 
(panel (b)) as a function of the cut-off ajf with the average density n = 1. 

of local maximum in the T c — U curve, denoted by U max , and the critical temperature at 
U = 0, denoted by T c0 , as a function of oim in Fig- 121 We can see that, as the cut-off olm 
increases, the positions of local maximum U max move to small U and approach steadily to 
and 1.1 after olm > 6 for the slave fermion and the slave boson, respectively. The critical 
temperatures T c0 shown in Fig. EJb) initially increase largely when we increase the an, but 
their dependencies on it become very small after olm > 6 in both slave particle approaches. 
All of these behaviors show that, in the case of low density (e.g. n = 1 here), the finite 
cut-off such as «Af = 6 is a good approximation in all range of U. 

We now make a comparison between the phase diagrams in Fig. Q by these two slave 
particle approaches. One can see that the T c of the slave fermion is always smaller than 



2l| . the local maximum in the 



that of the slave boson. As mentioned in Refs. [l7[ and 
T c — U curve is unreasonable. We can see from Fig. Efa) that, when making a cut-off such 
as olm = 5, the position of the local maximum in the slave fermion picture moves to zero, 
which means the local maximum disappears. However, it seems that the local maximum in 
the slave boson picture could not be eliminated by merely making a large cut-off. There 
is still a maximum located around U max = 1.1 even when a large cut-off such as cvjvf = 9 
is made. Moreover, the position of the maximum almost does not depend on the cut-off 
when oiu > 6. This unsatisfactory feature may come from the approximation of relaxing 
the constrain of eq. (Hi, which has less severe effect on the slave fermion than on the slave 
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boson 2l|. On the other hand, the critical temperature T c q at U = should be identical 
to the Bose-Einstein condensation temperature T* deal of ideal Bose gas, which in the case of 



three dimensions can be determined by |26j, |31[ 

2-7T 

rpideal 

c " ¥ 



n 



-| 2/3 

(35) 



.#3/2(0). 

with n being the average particle density and gi{x) = X^m=i £ mx l m% ■ We have T tdeal = 1.10 
when n = 1 [32]. As shown in Fig. I2fb), T c0 approaches 1.18 and 1.66 for the slave fermion 
and the slave boson, respectively. Therefore, in this special case of n = 1, T c o of the slave 
fermion is more appropriate for it is closer to the expected value 1.10. As shown in the 
phase diagrams, both slave particle approaches yield the critical interaction U c = 6 when 
olm — 2, which slightly deviates from the well-known mean field value U c = 5.83. This can 
be attributed to the mean field approximation of eqs. (32) and (33) and the finite cut-off of 



21( . However, a remedy to the mean field conditions is very difficult 13, |2l| and, 
when including more types of slave particles, the calculation of T c at very low temperature 
is very hard too for the divergence and the multi-solutions at U > 6 [21]. Hence, when 
olm_ > 2, we actually do not work out T c in the very low temperature region in this work 



All the critical temperatures calculated above are concentrated on the commensurate 
state with fixed integer density (i.e. n = 1), which is important only to the exact quantum 
phase transition at T = 0. We next turn to investigate the T c q away from the integer 
filling and compare them with the critical temperature of ideal Bose gas. In Fig. |3J the 
dependence of T c0 on average density n is plotted by both slave particle approaches, and 
the critical temperature of three-dimensional ideal Bose gas is plotted by eq. (|35jl . We can 
see that the slave fermion curve is always below the slave boson curve and closer to that 
of ideal Bose gas. On the other hand, the deviations between both slave particle curves 
and the ideal gas curve become larger as the density increases. As mentioned above, the 
multi-occupation of slave particles on one site is allowed in mean field approximation. In 
the high density region, the multi-occupation may occur more frequently for one should 
take more types of slave particles into account. This is why the deviations from the ideal 
gas become large in this region. However, in the slave fermion approach, multi-occupation 
of the same type of slave fermion is excluded, which makes this approximation less severe 
than in the slave boson approach and gives a curve closer to that of the ideal gas. One may 
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FIG. 3: The critical temperature T c q at U = as a function of average density n with the cut-off 
oim = 6. The dashed and dotted lines correspond to the slave boson and slave fermion, respectively. 
The solid line is the critical temperature of ideal Bose gas in three dimensions. 



notice that there is an unsatisfied feature in the slave fermion picture too. The slave fermion 
curve and the ideal Bose gas curve cross at two different values of the filling, which implies 
that the slave fermion can not give the correct function dependence of T c0 on n. In the 
end, we would pay some attention to the mean field nature of our theory. As we know, T c q 
is calculated on base of the mean field approximation and hence applicable to any spatial 
dimensions. This will lead to the conclusion that the Bose-Einstein condensation occurs in 
one and two-dimensional ideal Bose gas at finite temperature, which is obviously wrong for 



the violation of the Hohenberg theorem 



34j . However, this is a well known flaw: mean field 



theory usually breaks down in low dimensions due to the large quantum fluctuation 



IV. DENSITY AND COMPRESSIBILITY 



The central features of the Mott-insulator at T = are the integer filling factor n a = 5 a ^ a i 
and the zero compressibility k = dn/dfi = 0. At finite temperature, there is no exact 
Mott-insulating phase because the local filling factor may deviate from an integer and the 
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FIG. 4: The average particle density n as a function of p, with U = 11 at various temperatures, 
obtained by the slave fermion (Fig. E{a)) and slave boson (Fig. Etb)) approaches. The fl and U 
are the dimensionless chemical potential and on-site interaction, respectively. Six types of slave 
particles (from n° to n 5 ) are taken into account in the calculation. The compressibility can be read 
out from the slope of the curve. 

compressibility k does not equal zero. We will investigate the finite temperature properties 
of these parameters in this section. In principle, when given the value of U, one can calculate 
every type of occupation numbers n a as a function of fi by combining eq. (f!?4"|) with eq. 
Then the average particle density n, which is a function of Jx too, can be obtained according 
to eq. (jHSJ)- After depicting the n — Jjl curve, one can read out the compressibility from the 
slope of the curve. In Fig. 0J we show the calculated n — p, curves with U = 1 1 at different 
temperature in both slave fermion and slave boson pictures. As the figure shows, there 
are "steps" on the curves at very low temperature, where the density is very close to an 
integer and the compressibility is almost equal to zero. It is therefore reasonable to call these 
regions the " Mott-insulator" at finite temperature We can see that the Mott-insulating 
regions in both slave particle pictures dimmish and disappear gradually as the temperature 
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is increased. However, the quantitative behaviors of them are a little different, i.e., at the 
same temperature, the compressibility of the slave fermion deviates from zero more greatly 
than that of the slave boson. This implies that the influence of temperature is greater on the 
slave fermion than on the slave boson. In Fig. 01 there are only two "steps" corresponding 
to n = 1 and n = 2. One can obtain the higher "steps" by including more types of slave 
particles and working on the larger p. Note that, although the qualitative behavior of this 
"step" structure is quite in accord with the zero temperature density profile the value 
of the density deep in the superfluid region is unreliable for the invalidity of perturbation 
theory. 



V. SUPERFLUID DENSITY 



As mentioned above, the action near the critical point can be expended in powers of the 
order parameter $o,0; 

S E , eff = f3Q - G-\0, O)|$ ,o| 2 + a 4 |$o,o| 4 . (36) 

If the coefficient a 4 of fourth order term is positive, the superfluid density n can be deter- 
mined by 

-o = |($o,o)| 2 = ^^, (37) 



which minimizes the action. We calculate the coefficient a 4 following the steps in Ref. 
and show the detailed calculation in Appendix [X] The final result is 



17] 



2(a +!)(« + 2) ^ Q 2(q + !)(<* + 2) Q+2 



(jjl - aU) 2 [{2a + 1)U - 2/i] - {a + l)U} 2 [{2a + 1)U - 2fj] 

2(a + l)(a + 2)U _ 2(a + l)(a + 2)/3 ^ a+u ^ ^ 



-n a+l - i — ,K '—^ n a+1 (l =f n c 



(fi - aU) 2 [fi - (a + 1)U} 2 (ji - aU)\n - (a + l)U\ 

where =F correspond to the slave fermion or the slave boson and n a is defined in eq. (|24j) . 
In Fig. we plot the dimensionless a 4 = L(3a^/ ' zt as a function of the chemical potential p, 
with U = 12 and = 4. We can see that the sign difference in a 4 (i.e., =f) leads to a very 
striking result, i.e., a 4 of the slave boson is always negative. This means the Landau free 
energy is not minimized but maximized at eq. (|37|h that is, the mean field theory of the slave 
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FIG. 5: The dimensionless fourth-order coefficient 04 = Lfia^jzt as a function of chemical potential 
\x by the slave fermion and the slave boson approaches. In the calculation, U = 12, ztf3 = 4, and 
four types of slave particles (n°, n 1 , n 2 , and n 3 ) are taken into account. 



boson is not stable and can not give the Landau second order phase transition. The physical 
reason for this instability may come from the condensation of the slave bosons due to the 
relaxation of the constraint of one salve boson per site. To see this point more explicitly, we 
depict in Fig. [HI the dependence of the superfluid density | $0,0 1 2 on the chemical potential 
p, at various temperature. The inset shows the value of (7(0,0)/2a4 as a function of /2 in 
the slave boson picture. We can see that it is opposite to the standard mean field theory at 



T = |5|,ll6j], e.g., the "superfluid density" grows up at the places where are Mott-insulating 
regions in the zero temperature phase diagram |35j ] . 

In the slave fermion approach, 04 is positive definitely. Then the Landau second order 
phase transition theory may be safely applied. As shown in the Fig. El the superfluid regions 
are consistent with those in the zero temperature phase diagram, and become small and 
disappear gradually when the temperature is increased. It is notable that the perturbation 
theory used here is valid only near the critical point. Thus the value of the superfluid density 
far away from the critical point may be incorrect. The accurate value can be obtained by 



means of the Bogoliubov theory 



2fJ, but this is beyond the scope of this paper. Because 



only four types of slave fermions are considered, the superfluid phase in the larger \x region, 
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0.4 




FIG. 6: Superfluid density |<3?o,o| 2 as a function of chemical potential \x with U = 12 at different 
temperature. Four types of slave particles (n°, n , n 2 , and n 3 ) are taken into account. The solid, 
dotted, dashed-dotted-dotted, and dashed lines are obtained by the slave fermion approach and 
correspond to zt(3 = 7, 4, 2.5, and 1.8, respectively. Inset shows the value of G(0,0)/2a4 in the 
slave boson picture with zt(3 = 7, which should not be regarded as the superfluid density because 
the 04 is always negative. 



which corresponds to the higher filling factor, could not be obtained in our calculation. 



VI. THE EXCITATION SPECTRUM 



The excitation spectrum of the quasiparticle and quasihole can be determined by the pole 
of the Green's function, that is, by the equation G _1 (k, uj) = 0. From eq. (J2EJ), we have 

. n a+l -n a 1 
> a + l - = — . 39 

It is easy to show that both slave particle approaches give the same excitation spectrum 
at T = for n a = with a' being an integer filling factor. At nonzero temperature, 
one should take more than one type of slave particle into account and eq. (|H9j) will have 
more than two solutions In the case of low density n — 1, it is reasonable to take only 
three types of slave particles (n°, n l and n 2 ) into account, that is, only the processes in 
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which the occupation of a site changes among n°, n 1 and n 2 are considered (l7|. By this 
approximation, we can obtain two low-lying excitation spectra analytically, 



In FigEl we show the excitation energies (u+fi)/ zt as a function of U at k = for various 
temperatures. We can see that the tip of the lobe, where the gap for quasiparticle-quasihole 
excitation disappears, moves to the small U region as the temperature increases. This picture 
is qualitatively consistent with the superfluid- normal phase diagram obtained in Sect. IIII1 
Another feature of this figure is that the energy gap for quasiparticle-quasihole excitation 
is enlarged as the temperature or the interaction U is increased. Recently, Konabe et al. 
have obtained the same result by a standard basis operator method Q. When comparing 
panel (a) with panel (b), we can see that, at the same temperature, the lobe by the slave 
fermion evolves away from zero temperature lobe more greatly than that by the slave boson, 
from which we can conclude again that the influence of temperature on the slave fermion is 
larger than that on the slave boson. Note that, except the point of U = 0, the quasiparticle 
and quasihole branches do not meet at finite temperature, e.g., the dotted line in Panel (a) 
and dashed-dotted-dotted line in Panel (b) (the gap is very small). This may be due to the 
finite cut-off approximation and can be remedied by including more types of slave particles. 
In addition, the two branches of spectrum always meet at (0,0) at high temperature. The 
reason is that, when the temperature is high enough, one would have n° = n 1 = n 2 = 0.333 
at U — and then the right hand side of eq. (jjOj) would always equal zero. In all, eq. (j4*Uj) 
is valid only at very low temperature and near the Mott-insulator, where the approximation 
of ocm = 2 is justified. 

By using the slave fermion approach, we show, in Fig. |H1 the dispersion relation of a two- 
dimensional atomic gas in a square optical lattice at various temperatures. The behavior of 
the dispersion curves at zero temperature, marked by the solid lines, is qualitatively similar 
to that obtained by Sengupta et al. 9] using a strong-coupling expansion approach. At finite 
temperature, we can see that, in the different regions of the Brillouin zone, the influence of 
the temperature on the gap for quaiparticle-quasihole excitation is different. In the region 
between (tt, 0) and (7r/2,7r/2), the gap increases as the temperature is increased. However, 
in the region from (tt/2, 7r/2) to (it, 0), the gap diminishes when we increase the temperature. 



w k + A* = 2 + 2 ek ( 2 ™ 2 _ nl ~ n °) 



(40) 
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FIG. 7: The excitation energies (u + fi)/zt as a function of the dimensionless on-site interactions 
U at k = for various temperatures. Panel (a) and (b) are obtained by the slave fermion and 
the slave boson approaches, respectively. The inset in panel (a) shows the energy gap between 
quasiparticle and quasihole excitations. The average density is set to be n = N/L = 1 and three 
types of slave particles (n°, n 1 , and n 2 ) are taken into account. 

One can see that the temperature affects the excitation spectrum greatly at the points (0, 0) 
and (tt, 7r), but at other points such as (ir,0) and (n/2,n/2) the influence of temperature is 
very small. 
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FIG. 8: The dispersion relation (u> + n) / zt of a two-dimensional atomic gas in a square optical 
lattice at various temperatures. The results are obtained by the slave fermion approach and with 
U = 7.5. The solid, dotted, dashed-dotted-dotted, and dashed lines correspond to zt(3 = infinity, 
2, 1, and 0.5, respectively. 



VII. CONCLUSIONS 



The finite temperature properties of the Bose-Hubbard model were investigated by both 
slave boson and slave fermion approaches. Many physical quantities were calculated in the 
mean field level and we found in general there are no qualitatively differences either by 
using the slave boson or slave fermion approaches. However, when studying the stability of 
the mean field state, we found that in contrast to slave fermion approach, the slave boson 
mean field state is not stable. The mean field phase diagram of superfluid-normal transition 
was obtained by making a relative large a cut-off (e.g. ckm = 6), and the effect of the 
finite types of slave particles was estimated. The unreasonable local maximum in the phase 
diagram can be eliminated in the slave fermion approach by increasing the cut-off, but it 
can not be eliminated in the slave boson approach. In the low density region, the critical 
temperature at U = by our mean field approaches is quite close to the Bose-Einstein 
condensation temperature of ideal Bose gas in three dimensions. The particle density was 
derived and depicted as a function of chemical potential fl. The " Mott-insulating" phase 
at finite temperature, where the filling factor is very close to an integer, was identified and 
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how it evolves as the temperature increases was demonstrated. The superfluid density was 
calculated and plotted as a function of chemical potential /2. We showed that the slave 
boson approach could not give a correct superfluid density because the coefficient a 4 in the 
Landau free energy expansion is always negative. The low-lying excitation spectra were 
obtained analytically by taking three types of slave fermions (n°, n 1 and n 2 ) into account, 
and were plotted as a function of U and in the k space at different temperatures. It was 
shown that, in the different region of the Brillouin zone, the influence of the temperature on 
the quasiparticle-quasihole gap is different. 
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APPENDIX A: THE COEFFICIENT OF FOURTH ORDER TERM 



In this section we turn to calculate the coefficient a 4 of fourth order term by both slave 

n 

fermion and slave boson approaches, following the steps in Ref . |17[ . We are only interested 
in the $00 terms and can write the partition function of eq. (fTljl as 

D§* 00 D$ Q0 Da a Da a e~ s , 



S = iL(3X + e |$oo| 2 + Yl a*,knM Q/ V kn , 

a/3 k,n 



(Al) 



where a a ^ n and a aj kn are the Grassmann variables in the slave fermion picture and the 
ordinary complex numbers in slave boson picture, M is a matrix: 

B Xi M 
Bt X2 A 2 
B 2 X3 A 3 



\ 



J 
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with 



A n = ^^e $ 00 , 



L/3 

-Da = 7^ _e O < I ) 007 



Xa = -ia>„ + c(a), 

in which c(a) is defined in eq. (|15jl. We then integrate the slave particle fields out of the 
partition function and obtain the effective action S' e ^, 

-<?' 



D<$>* m D® m e °«t 



D$* 0J D$ 00 exp 



iL/5A + e |$ 00 | 2 T^ ln [ det ^ M ] ) | 

k,n / J 



5^ = ^L/5A + e |$ 00 | 2 T^ln[det/3M], 



(A2) 



k,n 



where =F correspond to slave fermion or slave boson. The determinant of f3M can be written 



as 



det/3M = [HPxc 



£/3 XaX, 



1$ 



oo 



a+l 



In the case of small $ o, hi [det [3M] can be expanded by using 

1 



ln(l - A|$oo| +#|$ 



oo | 



A|$ 00 r + |(-2^ + 4 J B)|$oo| 4 + O(|$ 00 | & ). 



Then the effective action up to fourth order in $ o can be given by 



S' eff = iL(3\Tj2 Lln 11^ 



XaXa+l 



1$ 



00 I 



(A3) 



^ 4 Vv^ 



E 



ME^ 1 ) + E E 



X.aX.a+1 



4(a + l)(/?+l) 

XaXa+lX/3X/3+l 



1$ 



00 1 ) 



a [a-/?|>2 

where the summation over k leads to the number L of lattice sites. The coefficient 02 of the 
second-order term |$oo| 2 can be read out directly from eq. (|A3|) . 

, e o ( a + l "> 



eo - eo' 



XaXa+l 



E(« +1 



-fi + aU 



(A4) 
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where ± corresponds to slave fermion or slave boson. We can see that both slave fermion 
and slave boson approaches reduce to the same result — G _1 (0, 0). The coefficient a 4 of the 
fourth-order term |$oo| 4 can be rewritten as 



CLa = ± — 

2 



EE 



(a + l) 2 



+ 



2(a + l)(a + 2) 



XaXctXa+lXa+l XaXa+lXa+lXa+2 



(A5) 



where ± corresponds to slave fermion or slave boson. Typically, we calculate the first term 



1 \ - 



(a + lf 



XaXaXot+lXa+1 



4e 



(a + l) 



(a + l) 



(3 ^ [iuj n - c(a)) [iuj n - c(a + 1)]' 



(A6) 



The frequency summations in this term can be performed by using 



Res 



/(*) 



± 1 



+ Res 



c(a) 



1 



eP z ± 1 



(A7) 



c(a+l) 



Res[F(z)] c , s denotes the residue of F(z) at the pole c(a), and for the m-th order pole, can 
be determined by 

1 d m ~ l 



R«F(*)U)= lim 



- x [{z-c{a)TF{z)). 



=(«) zZ'cXa) (m - 1)\ dz r 
By using this equation, the value of eq. (jA6|) can be evaluated, which is 



(A8) 



(a + l) 



XaXaXa+lXa+l 



(a + l) 2 P 
(fi - all) 2 

2(a + l) 2 



[n a (l+n a ) +n a+l (l+n a+1 ) 
n a ). 



in 



(A9) 



(all -/i) 3 

After calculating the second term in eq. (|A5|) in the same way, we can obtain the final 
expression of 0,4 as shown in eq. (|H8|) . In the Mott-insulating region at T = 0, we have 
n a = 5 a ^ with g being the average density. For the slave fermion, the sign in eq. (|38|) is 
minus and then 0,4 is reduced to 



CI4 



6o 4 f g 2 (g + 1) 2 

L(3\[{g-l)U-^ + {^-gUf 

(g+l)(g + 2) g(g-l) 



- gUf[2^ - (2g + 1)U] [fM-(g- l)U?[(?g - 3)U - 2/x] 

- 9{9 + 1)U ] (A10) 

[{g-l)U-^-gUfj' 1 ) 

which is exactly the same as the results in Ref. [3| and Ref. [3]. However, for the slave 
boson, the sign in eq. is plus and the result is quite different Furthermore, as we 
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have shown in Sec. 13 is always non-positive, which leads to the instability of the slave 
boson mean field state. 
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